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Abstract 

Electroproduction of light vector mesons is analyzed on the basis of handbag fac- 
torization. The required generalized parton distributions are constructed from the 
CTEQ6 parton distributions with the help of double distributions. The partonic 
subprocesses are calculated within the modified perturbative approach. The present 
work extends our previous analysis of the longitudinal cross section to the transverse 
one and other observables related to both the corresponding amplitudes. Our results 
are compared to recent experimental findings in detail. 

1 Introduction 

In a previous work [Ij we analyzed electroproduction of light vector mesons {V = p", (p, uj) at 
HERA kinematics within the handbag factorization scheme which is based on generalized 
parton distributions (GPDs) and hard partonic subprocesses. The latter are calculated 
within the modified perturbative approach [2] in which the quark transverse momenta are 
retained. The emission and reabsorption of partons from the proton is still treated in 
collinear approximation. In the kinematical region accessible to the HERA experiments 
that is characterized by Bjorken-x (xbj) of the order 10~^, it is not unjustified to restrict 
oneself to the gluonic subprocess Y'g ~^ ^9 and the associated gluonic GPD . In a recent 
paper [3] we extended that analysis to larger values of xbj ( ^ 0.2) such as accessible to the 
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Figure 1: Typical lowest order Feynman graphs for the 7*5' Vg (left) and 7*g Vq 
(right) subprocesses of meson electroproduction. 

HERMES experiment, but restricting ourselves to the analysis of the least model-dependent 
amplitude, namely the one for transitions from longitudinal polarized virtual photons to 
vector mesons polarized in the same manner, —>■ VlP- This analysis necessitates the 
inclusion of the quark subprocesses 7*g Vq (see Fig. [T]) and the associated GPDs for 
sea and valence quarks. 

Here, in this work we are going to complete the analysis of vector-meson electroproduc- 
tion by studying the amplitude for transversely polarized photons and mesons, 7^p VtP- 
The basic idea of modeling this amplitude has been already described in Ref. [T| for the glu- 
onic contribution. The extension to the quark contribution is straightforward. The crucial 
point at issue with the transverse amplitude is that the quark transverse momenta which 
are retained in the modified perturbative approach in order to suppress configurations with 
large transverse separations of the quark and antiquark forming the meson, regularize the 
infrared singularities occurring in the transverse subprocess amplitude in coUinear approx- 
imation [H|5] at the same time. These special configurations seem to be responsible for the 
excess of the leading-twist contribution to the longitudinal cross section over experiment. 
Indeed, taking into account the quark transverse momenta fair agreement between theory 
and experiment is achieved p]. We are going to apply this modified handbag approach 
to vector-meson electroproduction for energies, W, in the photon-proton center of mass 
system (c.m.s.) between about 5 and 170 GeV and photon virtualities, Q^, between about 
3 and 100 GeV^ while xbj is less than <0.2. We will compare our results in detail with 
recent data from HERA, COMPASS and HERMES. 

The GPDs H for quarks and gluons which dominate the process of interest for unpolar- 
ized protons at small xbj, are constructed from the CTEQ6 parton distribution functions 
(PDFs) [6] through double distributions [H [S] • Applying the same model to the GPDs H 
we are also going to estimate the so-called unnatural parity amplitudes and to study their 
implications on spin density matrix elements (SDME) and the double spin asymmetry All 
describing the correlation of the helicities of the beam and target particles. 

In Ref. |9] the electromagnetic form factors of the nucleon have been used to extract 
the zero-skewness GPDs H, H and E for valence quarks. The forward limit of -E^ab 
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analog of the PDFs, determined in Ref. [9], can be utilized for the construction of E^^^ 
non-zero skewness with the double distribution model. With these GPDs at disposal we 
will also estimate the SDMEs measurable with a transversely polarized target as well as 
the transverse spin asymmetry Ajjt- 

The plan of the paper is as follows: In Sect. [2] we will sketch the modified handbag 
approach. In Sect. [3] the double distribution model for the GPD H is described. The 
results on cross sections and SDMEs obtained within the handbag approach are presented 
in Sect. [Hand [5], respectively. Sect. [6] is devoted to an estimate of the role of the GPDs H 
and Sect. [7] to that of the GPDs E. We will conclude this paper by a summary (Sect. [H]). 



2 The handbag approach 

We are interested in vector-meson electropro duct ion in a kinematical region characterized 
by large and large W but small XBj( ^ 0.2) and small invariant momentum transfer —t. 
In the handbag approach, the amplitudes for the process 7*p — *• Vp which can be extracted 
from vector-meson electroproduction applying the one-photon approximation, factorize 
into partonic subprocesses (see Fig. [1]) and GPDs comprising the soft, non-perturbative 
QCD. At large the amplitude for — > ViP dominates and factorization has been 
shown to hold for it rigorously [TUIIII]. The amplitudes for other photon- meson transitions 
are suppressed by inverse powers of Q. Besides the longitudinal amplitude we will consider 
only the transverse one, VtP in this work which is the most important one of 

the suppressed amplitudes at small —t. Proton helicity flip is suppressed by v^^/2m (m 
being the proton's mass) and can be neglected in calculations of cross sections and SDMEs 
obtained with unpolarized protons. In Sect. [7] we will however estimate the size of the 
proton helicity-flip amplitudes explicitly. It will turn out that these amplitudes are indeed 
small. 

In the region of small xbj the dominant contributions are provided by the GPD H. To 
the proton helicity non-flip amplitude the GPDs contribute in the combination 

where the skewness ^ is related to xbj by 

^^-^[l+ml/Q^]. (2) 

2 - XBj 

Here, my denotes the mass of the vector meson. The GPD E can therefore safely be ignored 
since it is not expected that it is much larger than the GPD H, see Sect. [7] where we will 
take up this issue again. The GPD H only contributes to the transverse amplitude and 
can also be neglected in calculations of the cross sections. In Sect. [6] we will return to this 
problem and estimate the size of its contribution. Since there is no parton helicity flip in 
the partonic subprocesses to the accuracy we are calculating them, the parton helicity-flip 
GPDs [12] do not occur. 
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The contributions from H to the 7*p Vp amphtudes read {i = g, q, Xg = 0, x, 



q 



(^) = ^ E ^aC^ t E ^' t = 0) i/^(a;, t) . (3) 



The first sum runs over the quark flavors a and denotes the quark charges in units of 
the positron charge e. The non-zero flavor weight factors, Cy, read 



c; = -c; = c = = q = i. (4) 

The explicit helicities in (|3]) refer to the proton while fi is the helicity of the photon and 
meson and A that of the partons participating in the subprocess. Only the t dependence of 
the GPDs is taken into account in the amplitudes ([3]). That of the subprocess amplitudes 
Ti provides corrections of order t/Q^ which we neglect throughout this paper. In contrast 
to the subprocess amplitudes the t dependence of the GPDs is scaled by a soft parameter, 
actually by the slope of the diffraction peak. 

There is a minimal value of —t allowed in the process of interest 

train = -^m^ • (5) 

As other effects of order (e.g. the GPD E) the tmin effect is neglected. We note in passing 
that our helicities are light-cone ones which naturally occur in the handbag approach. The 
differences to the usual c.m.s. helicities are of order my/^/W'^ [I2] and can be ignored in 
the kinematical region of interest in this work. 

The full amphtude is given by a superposition of the gluon and quark contributions 

M"" = M""' + M'''' (6) 

and it is normalized such that the partial cross sections for 7*p — > Vp read (A is the usual 
Mandelstam function) 

dt Mc\Tr{\m K(\An -n2 0{+)+,0(+)+ 



which holds with regard to the above-mentioned simplifications. The cross sections inte- 
grated over t are denoted hj aL and ar- The full (unseparated) cross section for 7*p — > Vp 
is 

(T = (Tt + eaL , (8) 

in which e is the ratio of longitudinal to transverse photon fluxes. The power corrections 
of kinematical origin given in Eq. ([2]) and in the phase space factor ([7]) are taken into 
account by us. With the exception of these kinematical effects hadron masses are omitted 
otherwise. In the expression ([7]) for the cross section the symmetry relation 

-M^;.',-,. = -M^:',,. (9) 
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has been used which is an obvious consequence of the definition ([3]) and parity conservation. 
This symmetry relation coincides with the one that holds for natural parity exchanges; we 
therefore mark this amplitude by a superscript N. Since the contributions from H to the 
amplitudes obey the relation 

M%^,^_^^ = -M^;^,^^^ , (10) 

there are no interference terms between A^^* and A^^* in the observables for unpolarized 
electroproduction of vector mesons as can easily be shown. The relation (JTOl) is also ob- 
tained for the exchange of a particle with unnatural parity. In analogy to the contributions 
from H the amplitudes related to H are marked by a superscript U. The contribution 
\ is neglected in ([7]). 

Let us now turn to the discussion of the subprocess amplitudes. As is well-known, 
for the kinematics accessible to current experiments, the handbag amplitude evaluated 
in collinear approximation, overestimates the longitudinal cross section although with the 
tendency of approaching experiment with increasing [H [131 El- Oiig may wonder 
whether higher order perturbative QCD corrections to the subprocesses may cure that 
deficiency. However, this does not seem to be the case. NLO corrections [TSl [IS] are very 
large due to BFKL-type logarithms ~ Inl/^ and cancel to a large extent the LO term at 
low and low x^y A recent attempt [17] to resum higher orders with methods known 
from deep inelastic lepton-nucleon scattering [T8| seems to indicate that the sum of all 
higher order corrections to the LO term is not large. In view of this unsettled situation 
it seems to be reasonable to proceed along the lines advocated in Refs. [H [3] by simply 
using the LO result and add power corrections, especially since such corrections are anyway 
needed in order to account for the large transverse cross section ax- Once the higher order 
perturbative corrections are better understood within the modified perturbative approach 
one may add them to the LO results. As long as they are of reasonable magnitude there 
is no principal difficulty in this. It may merely be necessary to readjust the GPDs and the 
meson wave functions appropriately. 

As in our papers [T|, [3] we will model the required power corrections by employing 
the modified perturbative approach [2] in the calculation of the subprocesses. In this 
approach the transverse momenta of the quark and antiquark, (defined with respect 
to the meson's momentum), entering the meson are kept. In contrast to the situation at 
the mesonic vertex, the partons entering the subprocess, are viewed as being emitted and 
reabsorbed by the proton collinearly. This scenario is supported by the fact that the GPDs 
describe the full proton and their dependence therefore refiects the proton's charge 
radius ((k^)^/^ ~ 200 MeV) while the meson is generated through its compact valence 
Fock state with a r.m.s. k^ of about 500 MeV fi9[ [20] . Instead of a meson's distribution 
amplitude allowance is to be made for a meson light-cone wavefunction \E'yL(r, /c^) in 
the modified perturbative approach [TSl [2D|- Here r is the fraction of the light-cone plus 
component of the meson's momentum, q', the quark carries; the antiquark carries the 
fraction f = 1— r. Quark transverse momenta are accompanied by gluon radiation. In Ref . 
[2| the gluon radiation has been calculated in the form of a Sudakov factor exp[— ^(r, b, Q^)] 
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to next-to-leading-log approximation using resummation techniques and having recourse 
to the renormalization group. The quark-antiquark separation, b, in configuration space 
acts as an infrared cut-off parameter. Radiative gluons with wave lengths between the 
infrared cut-off and a lower limit (related to the hard scale Q^) yield suppression, softer 
gluons are part of the meson wavefunction while harder ones are an explicit part of the 
subprocess amplitude. Congruously, the factorization scale is given by the quark-antiquark 
separation, fip = 1/b, in the modified perturbative approach. For more details we refer to 
Refs. pig. 

Since the resummation of the logarithms involved in the Sudakov factor can only effi- 
ciently be performed in the impact parameter space [2] we have to Fourier transform the 
lowest-order subprocess amplitudes to that space and to multiply them with the Sudakov 
factor there. This leads to 

= J c^rd26^^^(r,-b) J•;,^^,(a;,e,r,Q^b)a,(/i«)exp[-5(r,b,Q2)]. (n) 

The two-dimensional Fourier transformation between the canonically conjugated b and 
spaces is defined by 

/>) = j^Jd'k^exp[-tk^.h]f{k^). (12) 

The renormalization scale fiR is taken to be the largest mass scale appearing in the hard 
scattering amplitude, i.e. /ir = max (r(5,fQ, 1/b). Since the bulk of the handbag contri- 
bution to the amplitudes is accumulated in regions where /i/j is smaller than 3 GeV we 
have to deal with three active fiavors. A value of 220 MeV for Aqcd is used in the Sudakov 
factor and in the evaluation of from the one-loop expression. 

The hard scattering kernels JF* or their Fourier transform JF* occurring in Eq. (ITT!) , 
are computed from the pertinent Feynman graphs, see Fig. [1] The result for the glu- 
onic subprocess is discussed in some detail in Ref. [1] and we refrain from repeating the 
lengthy expressions here. For quarks, on the other hand, the hard scattering kernel for 
longitudinally polarized photons and mesons reads 




•^o+,o+ + ^o-,o- — ~^F\ irr -p [Ts — Tu] , (13) 



where Nc denotes the number of colors and Cp = {N^ — 1)/ (2A^c) is the usual color factor. 
For convenience we only quote the sum over the quark helicities since this is what appears 
in ([3]). The denominators of the parton propagators read 

T = i T = ^ (14:) 

In both the quark and the gluon propagators we only retain kx_ in the denominators of the 
parton propagators where it plays a crucial role. Its square competes with terms oc r{f)Q'^ 
which become small in the end-point regions where either r or f tends to zero. 
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While for longitudinally polarized vector mesons the quark transverse momentum is 
only needed for the suppression of the leading-twist contribution (the ^ limit), it 
plays an even more important role in the case of transverse polarization. In the coUinear 
limit the spin wavefunction of the meson is Ty oc ^'^^^(il) where ey denotes the polariza- 
tion vector of the meson. It however leads to a vanishing contribution to the subprocess 
amplitude since the number of 7 matrices in the Dirac trace [f] that gives the hard scattering 
kernel, is odd (see Fig. [1] ) 0. If one allows for quark transverse momenta a second term in 
the covariant spin wavefunction appears 

^^vKu = -^«7y,7.}i^^ (15) 

for which the number of 7 matrices in the Dirac trace is even. In eq. f|T5|) My is a soft 
parameter of the order of the vector-meson mass, e.g. a typical constituent quark mass. The 
transverse momentum four- vector K = [0, 0, k^] is suitably defined as the quark-antiquark 
relative momentum and represents one unit of orbital angular momentum in a covariant 
manner [T|. Expanding now the hard scattering kernel for transversely polarized vector 
mesons, one finds 

^;,,+,(e,x,r,Q2,ir) = ^;,,+,(e,a;,r,Q2,ki) + A^|5;,+,(e,x,r,Q2,ki)ir, 

+ AJ-X,(e,x,r,Q2,ki)ir,i^^ + --- . (16) 

Higher order terms in Kfj_ are neglected. In the spirit of the modified perturbative approach, 
the k^ terms in the propagator denominators are kept as has been done for the longitudinal 
amplitude. As already mentioned the first term in Eq. (ITBI) . generated by ^'^y{±l), vanishes 
as a consequence of the number of 7 matrices in the trace and, evidently, the second one as 
well after integration on k^. Hence, the third term in f|T6|) is the leading one for transversely 
polarized mesons and leads to the subprocess amplitudes 

= -'-fJdr ^.^l^VTir, ki) A^r+A , (17) 

in which g± is the transverse metric tensor 0. Note that the wavefunctions for longitudinally 
and transversely polarized vector mesons are different in general. The transverse amplitude 
is of order ]<.]_/ (MyQ). Noting that {'k']_y^'^/My is of order unity, one realizes that the 
transverse amphtude is suppressed by 

A^+,,+,cx (ki)i/VQ, (18) 

For the quark subprocess the trace includes the hadronic matrix element that defines the GPDs and 
is oc 7+ or oc 7"'' 75. 

Note that for longitudinally polarized vector mesons e(0) ~ q'/my. Hence, the spin wavefunction is 
oc ^' in this case up to mass corrections. A mass term oc my^v in the spin wavefunction for transversely 
polarized vector mesons has been investigated [4] . In coUinear approximation this term leads to an infrared 
singular twist-3 contribution. Such mass terms are neglected by us. 
^AU its elements are zero except = = — 1. 
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with respect to the one for longitudinally polarized vector mesons. 

Working out the kernels, one finds after summation over the parton helicity that the 
kernel for transverse photon and meson polarization is obtained from the longitudinal one, 
Eq. (fT3ll . by the replacement 



T — T — '^J- *5 

■J- Q -J- 11 



2 My 

The new propagator denominators read 



(19) 



We note in passing that Tg and represent the denominators of the gluon propagators in 
the LO subprocess 7*? — > Vq while Tq and Tb belong to the quark propagators (see Fig. 
[1]). With the help of partial fractioning [i = s, u, j = a, b) 

TiTj 2 CiTi -\- CjTj , (21) 

we can cast the Fourier transform of the transverse subprocess amplitude into exactly the 
same form as for the longitudinal one, Eq. (fTTl) . The kernel is then a linear combination 
of four Fourier transformed propagators. 

The denominators of the parton propagators in f[T^ are either of the type 

or ^ 

^ "^^d^i^TOQ^ ' ^^^^ 

where di > 0. The Fourier transforms of these propagator terms can readily be obtained: 

Ztt * 
1 



T2 = —Ko{^d2{±^-x)bQ) e{±^~x) 
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+ -H^^^^^d,{x±ObQy{x±0, (24) 

where Kq and Hq^^ are the zeroth order modified Bessel function of second kind and Hankel 
function, respectively. 

3 The double distribution model 

As in Refs. [H [3] the GPDs are constructed from the PDFs with the help of double distri- 
butions [71 '8] . Since this construction is described in detail in our previous papers we will 
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only recapitulate a few basic elements of it. For details we refer to Refs. [H [3]. The main 
advantage of this construction is the warranted polynomiality of the resulting GPDs and 
the correct forward limit 0. It is well-known that, at low x, the parton distribution 

functions behave as powers 6i of x. These powers are assumed to be generated by Regge 
poles pn |22] . We generalize this behavior of the PDFs by assuming that the t dependence 
of the double distributions and hence the GPDs are also under control of Regge behavior. 
Linear Regge trajectories are assumed for small —t 

ai = ai{0) + a[t , i = g,sea.,val , (25) 

with 6i = ai{0) for quarks and, as a consequence of the familiar definition of the gluon GPD 
which reduces to xg{x) in the forward limit, 6g = ag(0) — 1 for gluons. The trajectories 
are accompanied by Regge residues assumed to have an exponential t dependence with 
parameters 6j. The following ansatz for the double distributions associated with the GPDs 
Hi is therefore employed (cf. [3j) 

f(3 a t) = e^»* I 3 hm + [jl - Wlf - 

MP,a,t) e \p\ ^^^^^ 22-^+1 r2(n, + l) (1 - ' ^ ^ 

where 

hgiP) = \P\gm) rig = 2, 

hteM = gsea(|/3|) Sign(/?) n^ea = 2 , 

him = gvai(/5)e(/3) n,ai = l. (27) 

For the decomposition of the double distribution into valence and sea contribution we 
follow the procedure proposed in Ref. [1^ and write 



fUP,<y,t) = nP,a,t)Q{P) - f%-(3,a,t)Q{-P)- (28) 

In the forward limit, — 0, this decomposition is conform to the usual definition of sea 
and valence quark PDFs. 

The GPDs are related to the double distributions by the integral 

/■I r^-m 

Hi(x,^,t) = d(3 da6(f3 + ^a-x)fi((3,a,t). (29) 

J-l J-l+\l3\ 

For convenience we employ an expansion of the PDFs {(3 > 0) 

= p-'^ (1 - P) ^ , (30) 

j=0 

which is particularly useful at low (3 and allows to perform the integral 0291) term by term 
analytically; its use is also convenient if the integral (!29|) is carried out numerically. The 
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factor (1 — Z?)^"'"'"^ serves for canceling the corresponding factor in Eq. (l26l) and has the 
additional welcome feature of roughly accounting for the (3-^1 behavior of the PDFs. 
The ansatz (|30l) results in a corresponding expansion of the GPDs 

3 

H,{x, e, t) = e"^' H,,{x, t) . (31) 

j=0 

The definition of the GPDs is completed by the relations 

Hai~x,^,t) = H^{x,i,t), HU-x,i,t) = -HUx,^,t), (32) 

and 

HUx,^,t) = -l<x<-e. (33) 

In Eq. f l2^ so-called D terms for the gluons and the flavor singlet quark combination are 
ignored [23] • The D terms ensure the appearance of the highest powers of the skewness in 
the moments of the GPDs. They only contribute to the less important real parts of the 
amplitudes since their support is the region — ^ < x < ^. The corresponding imaginary 
parts are related to the GPDs at a; = ^(1 + 2k^/(rQ^)) which lies outside the support of 
the D terms. We take this in vindication of neglecting the D terms. 

For the expansions of the PDFs we will use the same parameters as in Ref. [3] with 
the exception of a little change. ZEUS |23] now provides data on the cross section for p 
production up to = 100 GeV^ whereas in [3] the fits to the CTEQ6M PDFs were made 
for < 40 GeV^. As one may check, the CTEQ6M gluon and sea quark PDFs are not 
well described above 40 GeV^ by the expansion quoted in Ref. [3] but the addition of a 
L^-term to 6g 

5g = 0.10 + 0.06L - 0.0027 , (34) 

improves the fits to the gluon and sea quark PDFs considerably as can be seen from Fig. 
[2] (L = In QVQ2^ Qo = 4 GeV^). The term is irrelevant below 40 GeVl The valence 
quarks are not needed at high Q^. The parameters of the expansions fl30l) are quoted in 
Tab. [H We stress that with the exception of 6g, they are identical to those used in Ref. [3]. 
In the quoted ranges of Q2 and (3 the fits to the PDFs agree very well with the CTEQ6M 
solution; they are always well inside the band of Hessian errors quoted in [6]. Larger values 
of (3 are irrelevant to us since the region 0.5 < f3 affects the real parts of the amplitudes 
only marginally; the contributions are less than 0.5%. From the fitted PDFs the GPDs are 
evaluated with the help of Eq. (129!) . 

In an attempt to keep the GPD model simple we assume 

-^sea — -^sea — -f^s'ea 5 (35) 

where the flavor symmetry breaking factor is parameterized as 

Ks = 1 + 0.68/(1 + 0.52 In Q^Ql) , (36) 
as obtained from a fit to the CTEQ6M PDFs. 
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gluon 


strange 


•-^vai 


'-*'vai 


5 


Eq. (El) 


1 + 5o 


0.48 


0.48 


a' 


0.15 GeV"2 


0.15 GeV"2 


0.9 GeV"2 


0.9 GeV"^ 


Co 


2.23 + 0.362 L 


0.123 + 0.0003 L 


1.52 + 0.248 L 


0.76 + 0.248 L 


Cl 


5.43 - 7.00 L 


-0.327 - 0.004L 


2.88 - 0.940 L 


3.11 - 1.36L 


C2 


-34.0 + 22.5 L 


0.692 - 0.068 L 


-0.095 L 


-3.99 + 1.15 L 


C3 


40.6 - 21.6 L 


-0.486 + 0.038 L 









Table 1: The parameters appearing in the expansion (130|1 of the PDFs (L = In Q'^/Qq, 
Ql = A GeV^). 

As in Ref. [3] we take for the slope of the gluon trajectory the value a'g = 0.15 GeV~^. 
Since the sea quark PDFs are mainly driven by evolution for ^ 4 GeV^ it is assumed 
that asea(^) = ttg(t). A standard trajectory is adopted for the valence quark Regge pole - 
Q!vai(^) = 0.48 + 0.9 GeV~^ t. The parameter of the gluon residue is fixed from a fit against 
the HERA data for p [25] and 4> production [26]: 

2 

bg = 6sea = 2.58 GeV-2 + 0.25 GeV"' In ^ , (37) 

+ in 

The p and (p slopes of the cross sections practically fall together at HERA energies; there 
are only minor differences at low Q^. The parameter bg given in fl37|) leads to a t dependence 
of the differential cross section in perfect agreement with the recent ZEUS data |24j. The 
parameter of the valence quark residue is taken to be zero. This is in accord with the 
findings of the nucleon form factor analysis proposed in Ref. [9] in which the zero-skewness 
GPDs have been determined. 

We emphasize that the evolution of the GPDs is taken into account by us only approxi- 
mately through the evolution of the PDFs. This is reasonable since at low ^ the imaginary 
part of the gluon (and sea quark) contribution dominates which is oc H^{C,, ^, t) and there- 
fore approximately equals 2^g{2^) at low ^ (see, for instance [1]). Its real part as well as 
the valence quark contributions are only of importance near 4 GeV^, the initial value of 
the evolution. The very time-consuming numerical integration on a;, b and r forces us to 
use this approximative treatment of the evolution. In order to demonstrate the quality of 
our approximation we compare in Fig. [2] the gluon GPD at 40 GeV^ either obtained from 
our approximation or from evolving from the initial scale of 4 GeV^ using the evolution 
code developed by Vinnikov [27]. Only minor differences are to be noticed. 
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Figure 2: Left: The CTEQ6 gluon PDF at = 100 GeV^ as compared to our fit (solid 
line). The band indicates the Hessian errors of the PDFs. Right: The gluon GPD H'^ at 
^ = 0.01 and t = 0. Solid line: at the initial scale 4 GeV^. Dash-dotted (dashed) line: 
at the scale of 40 GeV^ evolved with the Vinnikov code [27] (approximate evolution). 



4 Results on cross sections 

Before we present our results obtained from the handbag approach we have to specify the 
meson wave functions used in the evaluation of the amplitudes. As in Refs. [H [3] we will 
take Gaussian wave functions {j = L, T) 



^vj{r,k^) = 87r2y2^/y,(/i^)a^^.[l + i?2''''(/i^)C'f'(2r-l)]exp[-a^^.kl/(rf)]. (38) 

Transverse momentum integration leads to the meson distribution amplitudes for which we 
allow for the second Gegenbauer moment besides the asymptotic form. The meson decay 
constants for longitudinally polarized vector mesons are known from the electronic decays 
of the vector mesons while those for transversely polarized mesons are taken from QCD 
sum rules [28]. In contrast to the decay constants fvL the latter ones are scale dependent 

fvTif^F) = fvTifio) [^^) ' ■ (39) 

Note that the decay constants for transversely polarized vector mesons always appear in 
the combination fyr/Mv, i.e. there is only one independent parameter in practice. In 
fact we use a typical constituent mass of 300 MeV for My and the QCD sum rule result 
[28] fyr/fvL — 0.8 at the scale of /iq = 1 GeV. The Gegenbauer coefficients have 
been found to be zero in the analysis of the longitudinal cross section [3j. Those for the 
transverse case are fitted to the data on the cxr or the cross section ratio R = olI^t- The 
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parameter 


PL 


Pt 




4>T 


fv[MeV] 


209 


167 


221 


177 


ay[GeV^^] 


0.75 


1.0 


0.70 


0.95 


-°2 


0.0 


0.10 


0.0 


0.10 



Table 2: The parameters appearing in the wavefunction fl38|) . quoted at the scale fio 
1 GeV. 



Gegenbauer coefficients are scale dependent 



,«s(/io)/ 



72j 



(40) 



where 72L = 50/81 and 72T = 40/81 [29]. Finally, the transverse size parameters ayj in 
fl38|l are either fitted to or to ctt depending on the polarization of the vector meson. 
The values for the various parameters are compiled in Tab. [2l Those for the longitudinal 
case are identical to the parameters used in [3j. 

The assessment of the theoretical uncertainties deserves special considerations. The 
results on cross sections (and other observables) are subject to parametric errors. The 
main uncertainties stem from the Hessian errors of the set of CTEQ6 PDFs. Since for the 
longitudinal cross section the parameters of the corresponding wave functions are adjusted 
such that good agreement between the data on ctl and the handbag results is achieved, 
there is no substantial additional uncertainties from the longitudinal wavefunction. Results 
for (Jl evaluated from sets of PDFs other than CTEQ6 also fall into the error bands in 
most cases (an exception is set for instance by the PDFs determined in Ref. [30l |3T] ) 
provided these PDFs are treated in analogy to the CTEQ6M set, i.e. they are fitted to the 
expansion fl30l) by forcing them to behave Regge-like with powers 5i as described above, 
and, if necessary, readjusting the transverse size parameters. The uncertainties in the ratio, 
i?, of the longitudinal and transverse cross sections mainly arise from the uncertainties 
of the wavefunctions for transversely polarized vector mesons (i.e. from the Gegenbauer 
coefficients, the transverse size parameters and from the ratio fvT/My). The errors due to 
those of the GPDs or PDFs which are mainly responsible for the theoretical uncertainties 
of the cross sections, cancel to a large extent in the ratio. 

The results for the longitudinal cross sections are the same as in Ref. [3]. We refrain from 
showing them again. The experimental data on the cross section ratio R are customarily 

by the relation 



determined from the measured SDME Tqq 



R = 



„04 
00 



e 1 



^04 
00 



(41) 



The SDME in fHTl) is understood to be integrated over the full range of t available in a 
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given experiment lI. The theoretical and experimental results on the ratio R are compared 
in Figs. [3] and |H In general we achieve very good agreement with experiment in particular 
with regard to the theoretical uncertainties displayed as shaded bands in the plots. 

The ratio R is mildly energy dependent for W larger than about 10 GeV while, for 
lower energies, it exhibits a somewhat stronger energy dependence, in particular at larger 
values of Q^, see Fig. HI This implies differences in the energy dependences of and ax 
which can be traced back to the different hard scattering kernels (see e.g. Eq. f|T9l) ) and 
the varying wavefunctions for longitudinally and transversely polarized vector mesons, in 
particular to the different values of the r.m.s. (~ l/ay)- An energy dependence of R 
seems to be indicated by the preliminary HERMES [35] and COMPASS |ii8j data although 
confirmation of this observation is demanded. Most of the older experiments have rather 
large errors such that a definite conclusion on a possible W dependence cannot be drawn at 
present. Unfortunately, for production there is only one low energy data point available 
P7] and this point is measured at the very low value of = 2.6 GeV^ which lies outside 
the range where the handbag approach, in its present form, can be trusted. We note that 
for (j) production the energy dependence of R is even milder than for the case of the p; the 
results for W = 10 GeV practically fall together with those at 75 GeV. 

In Fig. [5] the handbag results on the p and (p cross sections are compared to experiment. 
Again good agreement with the HI [25l[32] and ZEUS pHEHllM] data is to be observed in a 
large range of Q^. The leading-twist contribution to these cross sections, i.e. ctl evaluated in 
collinear approximation, is also shown. Although the leading-twist contribution approaches 
the experimental cross section with increasing there is still a small difference of about 
1.5(7 between both at = 100 GeV^ for p production. Note that even at that value of 
the transverse cross section which is included in the full one and also represents a power 
correction to the leading-twist result, is not negligible; it amounts to about 10%. The 
leading-twist contribution is about 20% larger than the one obtained within the modified 
perturbative approach at = 100 GeV^, i.e. the corrections due the quark transverse 
momentum have not yet disappeared completely. In Fig. [6] the energy dependence of the p 
cross section at a set of values is displayed. Within errors agreement is to be seen with 
the ZEUS data [J4J. This is not a surprise since the power 6g, related to the intercept of the 
gluonic Regge intercept, is fixed by the energy dependence of the HERA data on the cross 
sections, see Ref. [3]. As ai, see [5, the full cross section at HERA energies is dominated by 
the gluon contribution although the sea quarks are not negligible. Including the interference 
with the gluon the sea quarks contribute about 25% in the case of the and 40% in the 
case of the p at = 4 GeV^. The larger sea-quark contribution in the latter case is due 
to the flavor-symmetry breaking factor fl36p . Flavor symmetry breaking in the sea is 
important for the ratio of the and p cross sections, see Ref. [3] . Neglecting the sea quarks 
or assuming a flavor symmetric sea leads to an incorrect (f) — p ratio. Going to energies 
lower than about 10 GeV the valence quark contributions gradually become perceptible 
for p production. At, say, 5 GeV the valence quarks are responsible for about 40% of the 

^ The cross sections ctl and ctt have been obtained from integrating the differential cross sections over 
that range of t. 
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Figure 3: The ratio of longitudinal and transverse cross sections for p (left) and (right) 
production versus at W = 90 GeV and 75 GeV, respectively. Data taken from HI 
[25| [32l [33] (solid squares) and ZEUS [26l |3l] (open squares), respectively. The recent 
ZEUS data [23] are shown as open triangles. The solid lines represent the handbag results 
with the shaded bands indicating the uncertainties of the predictions. 




5 6 7 8 9 10 



Q'[GeV'] 




5 6 7 8 9 10 



Q'[GeV'] 



Figure 4: Handbag results for Ron p (left) and (right) production eitW = 75(10, 5) GeV 
shown as solid (dash-dotted, dashed) line. Data taken from Refs. |25| [26 | [3l]. Preliminary 
data from HERMES [351 [37] (sohd circle) and COMPASS j38] (diamond). The error bands 
are only shown at W = 5 GeV. For further notation refer to Fig. [31 
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cross section. We stress that the three contributions have to be added coherently. There 
are substantial interference terms which increase the cross sections markedly. For instance, 
at = 5 GeV and = 4 GeV^ the p cross section is doubled by the interference terms. 

A remark concerning the t dependence of the differential cross sections is in order. As we 
pointed out in [3] the cross sections drop exponentially with t to a very good approximation. 
Their slopes are approximately given by 2hg (see fl37|) ) plus a contribution from the gluonic 
Regge trajectory (see Eq. fl26l) and Ref. |3]). At large W the slopes of the differential cross 
sections for longitudinally and transversely polarized photons are the same while, at low 
W , the valence quarks generate small differences. Consequently, the ratio R is nearly t 
independent a fact that is in agreement with the ZEUS measurement [2^ . 

Finally, we check the theoretical consistency of the modified perturbative approach. 
Consistency is meant in the sense that the bulk of the perturbative contributions should 
be accumulated in regions where the strong coupling is sufficiently small. To find out 
whether or not this is the case, we set the integrand in fllip equal to zero in those regions 
where a<((/i/j) > ttg'^'*, and evaluate, say, the longitudinal cross section for p production as 
a function of af^^. Consulting Fig. [6] where the accumulation profile is shown, one observes 
that almost the entire result is accumulated in a comparatively narrow region of around 
0.4 at = 4 GeV^ and 0.25 at 40 GeV^. Thus the effective renormalization scales in these 
two cases are about 1.5 and 15 GeV^, respectively. Hence, our results on cross sections are 
theoretically self-consistent with regard to the above-mentioned criterion. Contributions 
from the end-point regions where the momentum fraction r tends either to zero or to one 
and where, in coUinear approximation, the renormalization scale becomes very small, are 
sufficiently suppressed. 



5 Results on spin density matrix elements 

In a number of experiments [25| [32| [38l [39] the SDMEs have been extracted from the 
measured decay angular distributions of the vector mesons. The formalism of the SDMEs, 
i.e. their relations to the amplitudes for the process 7*p — > Vp, has been developed by 
Schilling and Wolf [ID] long-time ago. This work has recently been repeated and extended 
to the case of a transversely polarized target proton by Diehl [41j. Since in the experimental 
papers the notation of Ref. is used throughout, we will adhere to it here in order to 
facilitate comparison. 

The SDMEs are given by properly normalized bilinear combinations of 7*]? — > Vp 
amplitudes. Due to the symmetry relations ([2]) and flTUl) there are no interference terms 
between M.^ and M.^ in the case of unpolarized electroproduction of vector mesons. The 
normalizations read 



2 

+ I ' 



Nr = 2 ^ [I A4t',++ r + I -^+.',++ 



(42) 



16 



4 6 



10 



20 
[GeV^] 



40 60 1 00 



10" 



10 

[GeV^] 



20 30 40 



Figure 5: The Q^- dependence of the cross sections for p (left) and (right) production at 
= 90 (p) and 75 GeV (p, 0). Data are taken from ZEUS ISSl [M] and HI [251 [32]. 
For p production at 1^ = 75 GeV the data and theoretical results are divided by 10 for the 
ease of legibility. The dashed lines represent the leading-twist contribution. For further 
notations refer to Fig. [3l 




Figure 6: Left: The cross section for p production versus W aX a. set of values. The data 
are taken from ZEUS [23] . For further notations refer to Fig. |3l The error band at 6 GeV^ 
is omitted for the ease of legibility. Right: The accumulation profile of the amplitude in 
dependence of a"'* at = 4 (solid line) and 40 GeV^ (dashed line) at ly = 75 GeV. 
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where helicity-flip 7* — > transitions are neglected as we do throughout this work. Up to 
a phase space factor, see Eq. ([7]), the normahzations are the differential cross sections for 
longitudinally and transversely polarized vector mesons. Since the SDMEs also provide a 
possibility to study the role of H we do not ignore the contributions from the corresponding 
amplitude M^yi in (H2il for later reference. 

Obviously, the SDMEs are functions of Q^,^ and W. Due to limitations in statistics 
it is not possible to measure the SDMEs as a functions of the three variables. Thus, 
frequently the SDMEs are presented as functions of one variable for average values of the 
other variables. For instance the SDMEs are quoted as a function of for an average value 
of W and integrated over all t available in a given experiment. Only in this case, and this 
is an important one to which we will mainly refer to in the following, the normalizations 
Nl and Nt refer to the respective integrated cross sections ai and cry up to the phase- 
space factor and eventually neglected suppressed amplitudes. If for instance the SDMEs 
are presented as a function of t the normalization represents differential cross sections at 
a fixed value of t but integrated over certain ranges of and W . 

The SDMEs allow for a separation of the absolute values of the two amplitudes, 0+ 
and A^++ ++, and for a test of their relative phases. Whether or not it is justified to neglect 
the other amplitudes can be examined as well. The SDMEs Tqq, r}_^ and lm.r'l_^ are related 
to the cross section ratio R. In terms of the amplitudes we are investigating these SDMEs 
read 

04 _ V I |2 



1-1 



-InirLi = E [I Ml,^,, P - I Ml,^,, r] . (43) 



Nt + eNl 

Neglecting the amplitude M.^ too these expressions simplify to 



= 2tU = -2lmrU = ■ (44) 



1 



eR 



In Figs. [3, El and [3 this prediction of our handbag approach is compared to experiment. The 
agreement is in general good within errors although with occasional exceptions. Thus, for 
p production, r|_^ and lmrl_i are not perfectly reproduced. On the other hand, their sum, 
measuring the double flip transitions 7!^ — >■ V-t, is nicely in agreement with zero. A further 
check is provided by the t dependence of the SDMEs. Since in the handbag approach the 
t dependence solely comes from the GPDs and these are identical for the two amplitudes 
the above SDMEs should be nearly flat in t which is indeed the case experimentally within 
errors [JUlSSj. Hence, the above three SDMEs are consistent with our assumptions and do 
not provide a significant signal for helicity-fiip amplitudes or contributions from H. The 
W dependence of the handbag predictions is displayed in Fig. [HI It is evidently very mild. 
In particular the results for the SDMEs at = 75 and 90 GeV fall practically together. 

After having checked that the absolute values of our two amplitudes are in fair agree- 
ment with experiment, we now turn to their relative phase Slt- In terms of our two 
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amplitudes the SDMEs relevant for the verification of the phase, read 



1 



lN kaN* 



2V2 1 + eR 



li 

cos 6 LT- (45) 



Predictions for the SDMEs in the case of production are shown in Figs. M and [TUl 
Fair agreement with experiments can be seen. For p production, on the other hand, a 
conflict is to be noted, see Figs. [7] and [HI The data [211 [251 [3S] require a rather large phase 
although with strong fluctuations (10 — 30°) while the handbag approach provides only a 
small value for it (e.g. 5lt = 3.1° at ly = 5 GeV and = 3 GeV^). Whether our model 
for the ^ amplitude which represents a power correction to the leading 7^ — >■ 
amplitude, is inadequate for this detail needs further investigation. However, that the sum 
Rer^o + Imr^Q amounts to only 1% of the corresponding difference of these SDME makes 
it clear that the neglected helicity flip '-/*—>■ V transitions in fH5l) are not responsible for 
the observed conflict. 

The HERMES collaboration has also measured the decay angular distribution of the p 
in the case of a longitudinally polarized lepton beam. These measurements which are not 
yet analyzed, will provide data on other SDMEs [IQ]. For instance, the SDMEs Imr^'g and 
Rer^Q also measure the phase 6lt introduced in (HSil 



ImrL = Rer?n = — ■= ^— sin^r^. (46) 

Predictions for these SDMEs are shown in Figs. [Til The other polarized SDMEs are only 
sensitive to the suppressed amplitudes, see Tab. [31 

In principle also the helicity flip '-/* V transitions can be calculated in the proposed 
handbag approach as well. But these transitions are strongly suppressed with respect to 
the amplitude A^o+,o+ as the amplitude for 7^^ — Vr transitions, see (fT8|) . In fact, 

7^ ^ 14 oc 

72 ^ Vr oc 



t 



Q 

t(ki)l/2 



2 



Q 

_f/U2 \l/2 

7t ^ V-T OC — , (47) 

where the powers of \/—i/Q follow from angular momentum conservation and the fac- 
tors (k^)^/^/(5 from the treatment of transversely polarized vector meson, see f[T8|) . The 
suppression factors given in P?]) engender the behavior of the SDMEs for helicity flip 
7* — > transitions listed in Tab. [31 The most important helicity-flip amplitude is the one 
for 7y Vl transitions. It is most clearly seen in Tqq as an interference term with the 
dominant longitudinal amplitude. It is definitely non-zero, of the order of 0.1 for p and 
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Figure 7: The SDMEs for the p meson versus 

(in GeV^) at W = 75 GeV. Data taken 
from [2U 1251 EH]- Note the different scales at the axis of ordinates. For further notations 
refer to Fig. [3l 
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Figure 8: SDMEs for the p meson versus 

(in GeV^) at W = 5, 10 and 75 GeV (dashed, 
dash-dotted and sohd hne, respectively). Preliminary data taken from HERMES [35] (solid 
circles) and COMPASS [38] (diamonds). Error bands only shown for W = 5 GeV. For 
further notations refer to Figs. [3] and HI 
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Figure 9: SDMEs for the d) meson versus 

(in GeV^) at = 75 GeV. Data taken from 
[26| [32] . For further notations refer to Fig. [31 




5 6 7 8 



0.3 



0.2 



0.1 



1-1 



Re 



5 6 7 8 




5 6 7 8 



Figure 10: SDMEs for the d) meson versus 

(in GeV^) at = 5 GeV. 
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Figure 11: The SDMEs Imr^Q and Rer^Q for p production versus at = 5 GeV. For 
further notations refer to Fig. [3l 



SDME 


amphtudes 


order 


'00; 'DO 

Re r 5*0, Rer^Q, Imr^Q, Imr^Q 

^00 

'"ii, ''"i-i, InirJ'_;^ 


1 (7t - V,) p 


-t/Q' 
v^(ki)/Q^ 

-t{K)IQ' 



Table 3: SDMEs controlled by helicity flip 7* — > V transitions. 

production in experiment [211 [251 [32] • At least HERMES [351 [26] observes a t dependence 
for it in agreement with expectation (oc v^^). The contribution of the 7^ — > ampli- 
tude to the cross sections, R and Tqq is at the percent level and can safely be neglected. 
The other SDMEs related to 7^ — >■ transitions are experimentally compatible with zero 
within errors. The only exception is to be seen in the recent high statistics ZEUS data [21] 
where Imr^Q and Rerj^Q have very small but non-zero for values of less than 10 GeV^. 
It would be interesting to see whether their t dependencies are in agreement with expecta- 
tion. For the remaining SDMEs being related to 72 — >■ Vr and 7!^ V-t transitions there 
is no significant deviation from zero at large experimentally. In view of these remarks 
it is fair to conclude that the neglect of the helicity flip 7* — >^ transitions in the cross 
sections for > 3 GeV^ is justified. 
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6 The role of H 

The expression for the amphtude is the same as that for M.^ given in ([3D except that 
the sum of the subprocess amphtudes is to be replaced by their difference 

■Kl,,^{V) = ^ E e,C^ dx [nll^^^ - nll^^_] H\x, t) . (48) 

The superposition of the various quark and gluon contributions is identical to that for the 
amplitude Ai'^ . The unnatural parity amplitudes satisfy the symmetry relation ffTOl) . It 
is evident from ( l48l) that parity conservation leads to a vanishing longitudinal amplitude 
-Mo+o+- The transverse subprocess amplitude is the same as in ( |T71) and ( fT9l) except that 
in the latter equation a plus sign occurs between TgTa and T„Tb now. 

The GPDs H are again modeled by the double distribution ansatz (126|) . With regard 
to the symmetry properties of H the functions hi now take the form 

hg{(3) = \(3\Agm) sign{(3) , 
hUP) = Ag,ea(|/3|), 

hUP) = AqU(3)Qi(3). (49) 

For the powers rii in the double distribution ansatz fl2^ the same values are taken as for 
the GPDs if*, see f l27l) . The decomposition of the double distribution into valence and sea 
contributions is made by 

/l(Aa,t) = [r{(3,a,t)-r{-f],a,t)]Qi(3), 

fUP, «, t) = f%P, «, t) 0(/5) + ri-P, «, t) e{-P) . (50) 

The double distribution ansatz for is incomplete because in moments of this GPD the 
highest power of ^ are lacking which leads to difficulties with the analytic properties of the 
amplitudes ^J]- We ignore this problem here since the contributions from seem to be 
unimportant. 

The required polarized parton distributions are taken from Ref. [43j and expanded 
according to 

km = r'- (1-/3) 2"-+^ E^.i/3^ (si) 

j=0 

using only integer powers. The resulting expansion parameters Cj and 6i are quoted in Tab. 
HI It is expected that the ai Regge trajectory controls the low-x behavior of the polarized 
valence quark PDFs. Since there are no recurrences of the ai(1260) we are forced to assume 
the standard value of 0.9 GeV~^ for the slope of the trajectory 0. Combining this with the 
spin of the ai, we obtain aa^(O) ~ —0.36 for the intercept. Such a low value is however 



^ Accepting exchange degeneracy for the ai and 772 (1617) trajectories the slope of the trajectory is fixed 
by the meson spectrum and is indeed 0.9 GeV^^. 
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vai 


val 


6 


0.48 


0.48 


0.48 


0.48 


Co 


0.61 + 0.033 L 


-0.320 - 0.040 L 


2.204 


-3.114 


Cl 


0.410 - 0.377L 


-1.427- 0.176 L 


-2.204 


8.096 


C2 


5.10 - 1.21L 


0.692 -0.068 L 


0.0 


-6.477 


C3 


0.0 


0.0 


0.0 


1.295 



Table 4: The parameters appearing in the expansion ( 130|) of the polarized PDFs and the 
forward hmits of E"-. The latter are taken from Ref. 0. The expansion (l30l) provides a 
fit to the Bliimlein-Bottcher PDFs [13] in the range 10"^ < P < 0.5 and 4 GeV^ < < 
40 GeV^. The powers S are kept fixed in the fits to the PDFs. 



in conflict with the small-x behavior of the polarized valence quark PDFs determined in 
Ref. 03], for which the power is rather about 0.7. As a compromise we therefore take the 
standard value of 0.48 for it. For the slope of the Regge trajectory we again take the value 
of 0.9 GeV^ and^for its residue 6vai = 0. 

The GPDs if* are obtained from the functions fi by the integral (1291) . They satisfy the 
relations 

H'{-x,^,t) = -H3{x,i,t), HlS-x,i,t) = HU^.^.t), (52) 

and 

Hl,{x,^,t) = -l<x<-e. (53) 

We checked that our proposed GPDs ifvai are in agreement with the data on the axial 
form factor for — 1<0.6 GeV^ and with the low —t (low x) behavior of H determined in 
[9]. In constrast to the situation for H, and H^^^ have opposite signs as a consequence 
of the behavior of the polarized PDFs. The lowest moments of the latter are known from 
P decays (see, for instance, Ref. [IS]). The usual assumption of a smooth behavior of 
the PDFs without a change of sign, leads to opposite signs of Au^ai and Arfvai- As our 
numerical studies reveal the gluon and sea quark contributions to Ai^ are very small and 
compensate each other to a large extent since the gluonic and sea quark GPDs H have 
opposite signs. Their combined contributions are practically negligible. This is the reason 
why we quote only the expansion parameters of hi for the valence quarks in Tab. HI 

Neglecting as in the preceding sections proton helicity flips and photon-meson transi- 
tions other than L —>■ L and T —>■ T, one can project out the unnatural parity amplitude 
for T ^ T transitions from a particular combination of SDMEs (see Eq. 



U 



1 r 



'00 ^'1-1 



^ Whether this is a consequence of lack of low- a; data in the PDF analysis or due to disregarded 
high-lying Regge cuts, is unknown at present. 
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Figure 12: Left: The ratio of au and a for p production versus Q"^ at W = 5 GeV. Data 
taken from HERMES The sohd (dashed, dash-dotted) hne represents our estimate 

(obtained with CuH^sx + CdH^^^y with CuH^g^ — edH'^^^), see text. Right: The hehcity cor- 
relation All for p production aX W = b (10) GeV dashed (dash-dotted) hne. Data taken 
from COMPASS and HERMES [19]. For other notations cf. to Figs. [3] and H 



This is the unnatural parity part of Nt scaled by Nt + eNl, see ( H2l) . Integrating over t one 
arrives at a cross section au defined in analogy to ctt in ([7]). Evaluating this cross section 
for p production from the amplitudes given in ([3]) and f HHl) and using the GPDs H described 
above, we find the results shown in Fig. [121 The cross section (Ju{p) is rather small but 
in agreement with the preliminary HERMES result [35] at = 2.88 GeV^ within an 
admittedly large error. For larger energies our approach will lead to even smaller values 
for au since the valence quark contribution disappears and, as we mentioned above, the 
combined gluon and sea contribution is very small (the typical size of gluon plus sea quark 
contribution to au{p) is 0.013 nb). We note that for p production, the HI data [25] provide 
values for au that are compatible with zero (e.g. at = 3 GeV^, au/c = 0.03 ± 0.07) 
while the ZEUS results |24j are about 1.5a above zero (e.g. at Q2 ~ 3.4 GeV^ au/a = 
0.03 ± 0.02). Both the experimental results are in agreement with our estimates within 
errors. An immediate consequence of the cancellation of gluon and sea contributions to the 
unnatural parity amplitude is that au for production is very small, in fact compatible 
with zero. This is in agreement with the preliminary HERMES data ^37j and with the HI 
data [32] . Thus, there is indication from both theory and experiment that au for p and 
production is small. Its neglect in ax seems to be justified. A larger cross section au is 
to be expected for u production because the combination euH^^^^ + ^dH^a,i occurs (see dl])) 
which is larger than e^iif^a^i — edH^^i given the relative sign of if^g^j and -ffyai- 

One may wonder whether or not it is possible to generate a value for au as large as. 
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say, the face value of the prehminary HERMES result [35] by using GPDs constructed 
from the polarized PDFs via the double distribution ansatz. In order to examine this 
issue we recall that the polarized PDFs are the difference of PDFs for helicity parallel and 
anti-parallel to that one of the proton while the unpolarized PDFs represent their sum. 
Suppose the gluon and sea quark contributions still cancel and assume now that the helicity- 
parallel distributions dominate which, in the limit x —>■ 1, follows from QCD [H]. In this 
case the double distribution ansatz leads to H^^^ = H^^^. Admittedly this is an extreme 
example since in all analyses [l3l HSl |16] the polarized d-quark distribution is negative. 
Another, more moderate example is set by the assumption H^^i = H^^i and H^^i = — i?vai> 
implying d quarks with dominantly antiparallel helicity. The results obtained from these 
two scenarios are also shown in Fig. [121 Obviously, it seems difficult to obtain agreement 
with the HERMES result [35] with GPDs constructed from the double distribution ansatz 
except extreme scenarios are realized in nature. 

The size of the amplitude Ai'^ can be elucidated further by considering the initial state 
helicity correlation A^l which can be measured with a longitudinally polarized beam and 
target. In contrast to cross sections and SDMEs where the corrections are bilinear in the 
terms and, hence, very small, the leading term in All is an interference between the 
Ai^ and the Ai'^ amplitudes. In fact, with the help of parity conservation as well as ([9]) 
and (fTOl) . one obtains (see pj) 



Re 



jvij — 1_ J — 1_ jvi^ — ^ _| — I 



AllIcp ^ eVp] = aVi^ ' ""'7 ""'""^ . (55) 



Nt + eNl 

We stress that here the target polarization is specified relative to the virtual photon direc- 
tion while in experiment it is usually defined with respect to the lepton beam direction. 
The conversion from our specification to the one used in experiments leads to a factor 
cos^^ in (!55|) [47j where 9^ denotes the angle describing the rotation in the lepton plane 
from the direction of the incoming lepton to the one of the virtual photon. This angle is 
given by |47j 



cos^, = 1 1-7^ , , \' ' ^l--f(l-y). (56) 




y - ^'7V4 ^ ^ _ 
1 + 7' 2 

The two parameters appearing in (156|1 are y = iW"^ + Q' — m^)/{s — m^), one of the 
conventional variables of electroproduction, and 7 = 2x-Q-^m/ Q. In the kinematical situation 
of interest 7 is very small and, hence, cos6'^ ~ 1. According to (ITS!) An is of order {k\) /Q"^ 
and, therefore, expected to be very small. In Fig. [12] our results for p production are shown 
at = 5 and 10 GeV. For the lower energy the valence quark contribution generates values 
for All of about 0.1 while at 10 GeV only extremely small values are found. The valence 
quark contribution has nearly disappeared at that energy and, as we mentioned above, the 
gluon and sea quark contributions cancel each other to a large extent. For instance, at 
W = h GeV and = 3 GeV', All{p) changes by —0.002 if the gluon and sea quarks 
are neglected. Because of the very small combined gluon and sea contributions we also 
predict very small values of All{.4')- Fo^' instance, at W = 5(10) GeV and = 3 GeV', 
Aii(0) = -0.002(-0.007). 
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Recently, the COMPASS collaboration has measured All for p production [H]. For 
g2 less than 2 GeV^ COMPASS observes very small values for All, which are compatible 
with zero. Their only data point for which its is sufficiently large for application of 
the handbag approach, is inconclusive because of its extremely large error; it is at variance 
with our results by about la, see Fig. [121 The HERMES results on this observable for 
p and (p production [19] is in agreement with our predictions. The SMC experiment [50] 
observes a double spin asymmetry for p production a.tW = 15 GeV that is compatible with 
our results within large errors. In passing we remark that All is sensitive to the relative 
phase 6nu between the amplitudes and A1++_,_+. Therefore, a large value of au 

is not necessarily in contradiction with a small value of All provided the phase is near 
90°. Nevertheless, in our approach the phase is small. For instance, at = 5 GeV and 

= 3 GeV^ we find S^u = 3.7°, i.e. a large value of au would go along with a large value 
of All in our approach. For instance, a scenario with H^^y = H^^y leads to All{p) = 0.14 
for this kinematical point of reference. 

Let us return to the issue of the size of the combined gluon and sea contribution to 
the unnatural parity amplitude. Up to know we assumed that this combined contribution 
is very small as follows from the double distribution ansatz using current polarized PDFs. 
One may wonder whether these PDFs are really correct or whether the smallness of the 
combined gluon and sea contribution is perhaps a special feature of our double distribution 
ansatz for H. First we note that in all current analyses of the polarized PDFs [^ Il5l H6] 
Ag and the polarized sea have opposite signs and are rather small in magnitude [j. In 
particular a large positive Ag is in conffict with measurements of the All asymmetry in 
the production of jets [51] and vr^ mesons [52] in inclusive proton-proton collisions. A 
negative polarized sea is for instance seen in the HERMES semi-inclusive deep-inelastic 
scattering data ^53j- Thus, we think that the main features of the polarized PDFs are 
correct. A small combined gluon and sea contribution from H seems to be required by the 
relevant data too. Leaving aside the HERMES results on au and All for p production 
which receive contributions from Hy^i, we note that all other pertinent data are small and 
in most cases compatible with zero. These data are au{p) from HI [25j and ZEUS [24j, the 
same cross section for (p production from HERMES [37j and HI |i32j and finally the All 
data from Compass [48j and SMC [50j. Thus, scenarios in which the combined gluon and 
sea quark contribution is large in magnitude seem to be excluded by the current data. Our 
double distribution ansatz for and Hsea , as naive it may be, qualitatively reproduces 
the main features of the data. 



7 Proton helicity flip 

The analysis of observables for vector-meson electroproduction measured with a trans- 
versely polarized proton target requires the proton helicity flip amplitude. The handbag 

^ A negative is also discussed in Ref. [45] . 
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contribution to this amplitude is given by 



K-,,+ iV) = -| ^ E eaC^ dx [nll^^^ + niL,,:] E\x, e, t) , (57) 

where tmin is ignored. Our choice of the phase of this amphtude is in accord with con- 
ventions exploited in ^T], H?]. In general there is also a contribution from the GPD E 
feeding the amplitude A4^!_ It is expected to be small and neglected in our estimate 
of the size of the proton helicity flip amplitude. The evaluation of Ai^ for proton helicity 
flip is analogous to that of the non-flip amplitude ([3]) except that the convolution is now 
to be performed with the GPD E instead of H. The construction of E through double 
distributions is also analogous to that of H, see Sect. El The only but crucial difference 
is that the forward limit e{x) = E{x, 0, 0) is inaccessible in deep inelastic lepton-nucleon 
scattering. However, the forward limits of the valence quark GPDs have been determined 
phenomenologically in the form factor analysis performed in Ref. [9]. The parameters of 
e"g^i and e^^^j expanded according to (ISTI) . are taken from [9]. They are quoted in Tab. H] 
at a scale of 4 GeV^. For other scales these functions are unknown which does not matter 
here since we will estimate proton helicity flip only for photon virtualities near that value. 
Note that in contrast to Uvai and d^ai, e.^ax and e^i have opposite signs. This is due to the 
fact that they are normalized by 

1 

cixe^^i(x) = Ka, (58) 







where k°- gives the contribution of quark flavor a to the anomalous magnetic moment of 
the proton («;„ ~ 1.67, ~ —2.03). The forward limits of E for gluons and sea quarks 
are unknown. But there is a sum rule 



/ dxxeHx) = -V / dxx el.Jx) - 2Y^ [ dxxe^ix), (59) 
Jo ^ Jo V ^0 

which follows from a combination of Ji's sum rule and the momentum sum rule of deep 
ineleastic lepton-nucleon scattering [12]. Neglecting a possible difference between and 
e^, we can evaluate the valence-quark term in the sum rule fl59l) from the GPDs specified 
in Tab. HI We obtain 

V / rfxxe^^i(x) = 0.008 ±0.007. (60) 
a ^0 

This signals a remarkable compensation between the second moments of e"^i and e'^^^ 



^ ~ 0.026, (61) 



^valJ 



which is even stronger than that of their first moments ((k„ + K,d)/{i^u — i^d) — 0.1). The 
error in ([HDD has been estimated from those quoted in P]. Hence, the moment of in 
(!60|) is only about as large as the sum of the sea quark moments. This is to be contrasted 
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with the situation for H where the corresponding gluon moment is more than four times 
larger than the sum of the sea quark ones. Another argument that points into the same 
direction is the behavior of the gluon (or Pomeron) Regge trajectory. As is well-known 
this trajectory couples mainly to the proton helicity non-flip vertex while the flip coupling 
is very small. It is hard to find phenomeno logical evidence for a non-zero flip coupling 
[5l] . Thus, the relative importance of gluon and valence quark GPDs is very different for 
E and H. It seems unlikely that plays an analogously prominent role as H^. In a flrst 
step we therefore assume that, for HERMES kinematics, proton helicity flip is dominated 
by the valence quarks (see also the discussion in Ref. [16]). Since the same Regge poles 
contribute to E and H, we therefore use the standard valence-quark trajectory here as 
well and assume bl^^y = 0, too. With regard to this situation we cannot estimate the size 
of proton helicity flip for production but we expect it to be very small. We stress that 
due to the opposite signs of E^^^ and E^^^, their contribution to p production off protons, 
oc CuEl^^i — CdE^g^y, is much smaller than that from the corresponding contribution of i^vai- 
This provides additional justiflcation for the neglect of E in the proton helicity non-flip 
amplitude (see discussion after Eq. ([2])). An interesting case is u production since for it 
the combinations CuE'^^^ + eaE^^^^ and e^H^gx + ^dH^^i occur. The flrst combination is larger, 
the second one smaller than for p production and, hence, a markedly larger ratio of proton 
helicity flip and non-flip is expected for cu production. For instance, at = 4 GeV and 
W = 5 GeV the flip/non-flip ratio of the absolute values of the uj amplitudes is about 13 
times larger than the corresponding ratio for p production. 

Recently the formalism for the SDMEs in the case of a proton target polarized per- 
pendicular ('normal') with respect to the plane in which the scattering 7*p Vp takes 
place, has been developed [41j. These SDMEs are denoted by n^^/ and related to bilinear 
combinations of the amplitudes for helicities p, p' and cr, a' of the virtual photon and the 
meson, respectively: 
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(62) 



If one neglects the amplitudes Ai^ in accord with our flndings described in Sect. [6l as well 
as helicity flip transitions 7* — y only the SDMEs 
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are non-zero. Explicitely these SDME read 
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For non-zero SDMEs n^^^, phase differences between the proton hehcity flip and non-flip 
amplitudes are mandatory. Such phase differences are provided by the handbag approach 
since the non-flip amplitudes are built up by gluons and quarks while the flip amplitudes 
receive only contributions from the valence quarks in our model for the GPD E. Indeed 
we obtain the values 38.8° and 34.7° for the phase between the proton flip and non-flip 
amplitudes in the case of longitudinal and transverse photon polarization, respectively. In 
Fig. [H the SDMEs ^ are shown versus at IV = 5 GeV, the trivial factor ^^jlm, 
see ( 1571) . is pulled out and t set to zero otherwise. The scaled SDMEs are small, of the 
order of five percent. These SDMEs will be measured by HERMES. 

One may also consider transverse proton polarization lying in the 7*p — > V'p plane 
('sideways'). In this case SDMEs, denoted by [H], occur that are analogous to (!62|) 

but with a plus sign between the two terms. The SDMEs s^^, for photon-meson helicity 
non-flip are given by products of two small amplitudes, M. for proton helicity flip and 
MP . They are therefore very small in our approach. 

Finally, we estimate the asymmetry Aut- of ep eVp for a transversely polarized 
target, normal to the 7*p Vp scattering plane. It is measured as the sin (0 — 0^) 
moment of the electropro duct ion cross section where (p is the azimuthal angle between the 
lepton and hadron plane and (ps the azimuthal angle of the target spin vector defined with 
respect to the direction of the virtual photon [51]. As for the asymmetry All the conversion 
of this spin vector into the one used in the experimental setup where the target polarization 
is specified relative to the lepton beam, again leads to a factor cos 9^ in principle. According 
to the discussion following Eq. (1561) it is replaced by one. In the handbag approach the 
dominant contribution to this asymmetry reads 

AuAep eVp) = 4 -— — (65) 

i\t + ^jVl 

It is just the imaginary part of the sum of and en^Q and is also proportional to y/^/2m. 
We again pull out the latter factor and display the scaled asymmetry, evaluated at t = 0, 
in Fig. [121 We obtain a positive asymmetry. In contrast to All it is finite to leading- 
twist order which is obtained from (165|) by neglecting the contributions from transverse 
photons and evaluating those from longitudinal photons in collinear approximation. For 
comparison the leading-twist contribution is also shown in Fig. [131 It is not too different 
from the full result. A preliminary HERMES result ^55j for p production, integrated on 
the range > -t > 0.4 GeV^ is -0.033 ±0.058 at = 3.07 GeV^ and 1^ = 5 GeV while 
we find 0.02 ±0.01 for this kinematical situation. Note that the scaled asymmetry is still t 
dependent although mildly. Ignoring this and integrating just one makes an error of 
about 10% at ^ 3 — 4 GeV^. For u production Aut is about 10 times larger than for p 
production. For (p production, on the other hand, we expect a very small asymmetry since 
the gluon and sea contributions are not only small but cancel each other to some extent, 
see dSS]). 

The asymmetry Aul for an unpolarized beam and a longitudinally polarized target is 
given by the same expression as Aut- Only the mentioned conversion factor cos^^ is to 
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Figure 13: Left: The imaginary parts of the SDMEs Uqq (sohd), n^'^ (dashed) and no+ 
(dash-dotted hne), scaled by 2m/v^^, for p production versus at W = 5 GeV. The 
error band is only shown for ngg. Right: The asymmetry Ajjt, scaled by 2m /y/^, for 
p production versus Q"^ at W = 5 GeV. The dashed line represents the leading-twist 
contribution. 

be replaced by sin^-y which is very small [17j. The beam asymmetry Am obtained with 
a logitudinally polarized beam and an unpolarized target is zero given that helicity flip 
7* — > transitions can be neglected. 



8 Summary 

Together with Ref. [3j this work gives an exhaustive description of light vector-meson 
electroproduction within the handbag approach for a wide range of kinematics reaching 
from the HERMES up to the HERA kinematical settings. Our handbag approach includes 
power corrections which suppress the leading-twist amplitude for 72 — > Vl transitions 
and allows for a calculation of the transverse amplitude 'y^ Vr- In order to specify 
our handbag approach fully we have to mention the soft physics input, namely the GPDs 
that are constructed from PDFs with the help of double distributions, and the light-cone 
wavefunctions for the mesons. GPDs and wavefunctions affect the handbag amplitudes 
differently and can therefore be disentangled. The wavefunctions provide effects of order 
(k5^)/Q^ controlled by the transverse size parameter ay while the GPDs mainly influence 
the XBj dependence of the amplitudes or, at fixed Q^, the W dependence. Besides the 
dominant contributions from the GPD H ('natural parity' contribution) we also estimated 
effects from the GPDs H ('unnatural parity' contribution) and E controlling the proton 
helicity flip amplitudes. These effects are generally small. With our analysis we achieve 
a very good description of the HERA, HERMES and COMPASS data on the separated 
and unseparated cross sections for p and electroproduction, on the ratio ctl/ctt, on the 
SDMEs and on some spin asymmetries. The only problem we observed is that the relative 
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phase between the longitudinal and transverse amplitudes seems to be larger in experiment, 
in particular in the HERMES experiment [35], than our handbag approach predicts. The 
neglect of contributions from transitions other than 72 — > Vl and — *■ Vr to the cross 
sections seems to be justified. Only little contributions from 7^ — *■ transitions are to be 
observed in some of the SDMEs experimentally. We note that in Ref. [56] the longitudinal 
amplitude has also been analyzed within the handbag approach. The main difference to 
our work is that in [56] the gluonic contribution is treated in leading-log approximation 
[57] and added incoherently to the quark amplitudes. This line of action understimates 
the gluonic contribution at low energies. 

The applicabihty of our approach is limited to > 3 — 4 GeV^, > 4 — 5 GeV and 
xbj ^0.2. The restriction of is due to the mentioned - and still unsettled - difficulties 
with higher-order perturbative corrections as well as due to the neglected corrections of 
order m? /Q"^ and —tjQ^. There may be also power corrections of other dynamical origin 
which become large at low . The restriction of W has its origin in the asymmetric 
minimum the cross sections exhibit at ~ 3 — 4 GeV. While the cross sections [2^ [251 
|3H [58] mildly increase towards larger W , they [35], [591 EO] increase sharply in the opposite 
direction. In fact they rise by nearly an order of magnitude between ~ 4 and 2 GeV. 
Obviously, a new dynamical mechanism seems to set in and the handbag physics is perhaps 
not applicable here. On the other side, it dominates for W >4 GeV. The mild increase 
of the cross sections with energy beyond the minimum is well described by the handbag 
physics as the results presented in this article and in Ref. [3] reveal. The restriction of xbj 
is of difference quality. It allows to neglect contributions of order x|j or e.g. in ([1]), 
which simplifies the analysis of vector-meson electroproduction strongly. 
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